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Abstract 

We study Bromine and Chlorine chemisorption on a Ag(lOO) surface, using a lattice-gas model 
and the quantum-mechanical Density Functional Theory (DFT) method. In this model the Br and 
CI ions adsorb at the fourfold hollow sites of the Ag(lOO) surface, which can be represented by a 
square lattice of adsorption sites. Five different coverages were used for each kind of adsorbate. 
For each adsorbate and coverage, we obtained the minimum-energy configuration, its energy, and 
its charge distribution. From these data we calculated dipole moments, lateral interaction energies, 
and binding energies. Our results show that for Br the lateral interactions obtained by fitting to 
the adsorption energies obtained from the DFT calculation are consistent with long-range dipole- 
dipole lateral interactions obtained using the dipole moments calculated from the DFT charge 
distribution. For CI we found that, while the long-range dipole-dipole lateral interactions are 
important, short-range attractive interactions are also present. Our results are overall consistent 
with parameter estimates previously obtained by fitting room-temperature Monte Carlo simulations 
to electrochemical adsorption isotherms [I. Abou Hamad et al., J. Electroanal. Chem. 554 (2003), 
211; Electrochim. Acta 50 (2005), 5518]. 



1 



I. INTRODUCTION 



The adsorption of halides on noble metals provides important model systems for studying 
adsorption on metal surfaces, particularly when there are ordered adsorbate structures. For 
this reason, these adsorption processes have been extensively studied |l, 2|. Adsorption of 
Bromine and Chlorine on metal has been the subject of many studies over the years. The 
systems we study here are Br and CI chemisorbed on single-crystal Ag(lOO). Experimentally, 
Kleinherbers et al. [3ij have found that the adsorption of Br, CI, and I on Ag(lOO) surfaces 
in vacuum all resulted in the formation of a c(2 x 2) overlayer with the adsorbates in the 
fourfold hollow sites. This implies a very strong, short-range repulsion, which we model as 
a nearest-neighbor exclusion [4]. 

The bonding of the adsorbates to the substrate and the surface electronic structures have 
been studied by Density Functional Theory (DFT) calculations. It is found that the bond 
between Br or CI and the substrate is covalent with a polarization due to electron transfer 
from the substrate to the adsorbate 5|-[7|. The polarization results in dipole moments on 
the surface, which cause long-range dipole-dipole interactions between the adatoms. 

Long-range dipole-dipole interactions have previously been incorporated in a lattice-gas 
model employed in room-temperature Monte Carlo simulation studies of the adsorbed sys- 
tem . In these works, the lateral interactions were extracted by fitting the results of 
the simulations to electrochemical adsorption isotherms. In the present study we instead 
estimate the lateral interactions by fitting the lattice-gas model to our DFT results. 

We extract the next-nearest-neighbor lateral energy and the binding energy by fitting the 
lattice-gas model to the adsorption energies obtained from the DFT calculation. The same 
DFT calculation also yields charge distributions from which dipole-dipole interactions can 
be directly calculated. By comparing the two results, we examine the significance of the 
long-range dipole-dipole interactions within the lattice-gas model. 

In this study we present DFT calculations using supercell models for the Ag(lOO) surfaces. 
The adsorbates in these DFT calculations are assumed to occupy a lattice of adsorption sites 
in accordance with a lattice-gas approximation 10|]. The lattice-gas assumption of strongly 
located adsorbates is consistent with previous DFT calculations and dynamic Langevin- 
equation simulations for a continuum model 

The adsorption energies and charge distributions were calculated by DFT. We assume 
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long-range dipole-dipole interactions between the adsorbates, and we implement these long- 
range interactions in the fitting of the lattice-gas model to the adsorption-energy results and 
the dipole moments obtained from the DFT calculations. 



'T ca^ 
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culations of adsorption 



13| and heteroepitaxy 



Estimates of short-range lattice-gas interactions from U. 
energies have also recently been performed for homoepitaxy 
systems. However, these studies do not consider charge transfer and long-range interactions. 

The rest of this paper is organized as follows. Section 2 describes the details of the DFT 
calculations and the methods used to analyze the results, Section 3 discusses the calculation 
of the dipole moment. Section 4 discusses the lattice-gas model. Section 5 presents the 
lattice-gas fitting, and Section 6 contains a discussion. 



II. DENSITY FUNCTIONAL THEORY 



We applied DFT to obtain the ground-state energies and electron density functions for the 
adsorption of Br or CI on a slab representing a Ag(lOO) surface. We prepared slabs with 
seven metal layers. Convergence checks with respect to the number of layers are discussed in 
Appendix A. The slab was placed inside a supercell with periodic boundary conditions. Two 
different sizes of supercells were examined. A 2 x 2 supercell with the size of 2ax2ax37.53A 
, and a 3 X 3 supercell with the size of 3a x 3a x 37.53A. Here, a = a/V2 where a = 4.17A is 
the lattice constant of bulk Ag, which we obtained from DFT calculations by minimization 
of an Ag fee structure. The 2x2 supercell contained four surface Ag atoms on each side of 
the slab (28 Ag atoms in total), while the 3x3 supercell contained nine surface Ag atoms 
on each side (63 Ag atoms in total). 

The orientation of the surface normal defines the z direction. To maximize the symmetry, 
we distributed the adsorbates on both sides of the slab. One, two, and three Bromine or 
Chlorine atoms were placed on each 3x3 surface to represent coverages 6 = 1/9, 2/9, and 
1/3, respectively. Two Bromine or Chlorine atoms were placed on each 2x2 surface to 
represent = 1/2 and one to represent 6 = 1/4. Here the coverage 6 is defined as 

where Q = 1 when the site is occupied by the adsorbate, and Cj = otherwise. In other 
words, the coverage is the number of adsorbates divided by the total number of all possible 
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FIG. 1: (A) The cross section of the superceh, (B) a three-dimensional representation of the 
superceh, and (C) surface distributions of the adsorbates for various coverages. Adsorbate atoms: 
gray. Surface Ag atoms: white. Bulk Ag atoms: black. 



adsorption sites, A^site- Figure [H shows the cross section of a supercell and surface distri- 
butions of the adsorbate for various coverages. Due to the nearest-neighbor exclusion and 
the periodic boundary conditions the adsorbates can only be placed in diagonal positions, 
limiting 6 to less than or equal to 1/2. 

The calculations were performed by the DFT method using the Vienna Ab Initio Simula- 
tion Package (VASP) 14-3 • 

The basis set was plane-wave, with the generalized gradient- 
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corrected exchange-correlation functional [l7|, [l8|, Vanderbilt pseudopotentials [l9|, |20|, and 



a cut-off energy of 400 eV. The fc-point mesh was generated using the Monkhorst method 
with a 7 X 7 X 1 grid for the 2x2 supercells and a 5 x 5 x 1 grid for the 3x3 supercells. 
To get to the configuration with minimum energy, we used a selective dynamics method, by 
which the ions in the top and bottom layers were allowed to relax in the z direction only, as 
opposed to the full dynamics in which the atoms would be allowed to move in all directions. 
This is the first step to avoid surface reconstruction, which is not expected to occur in this 
system under electrochemical conditions. The second step is to average the z coordinates 
of the top and bottom layers. The DFT results yield total energies and electron densities, 

We next ran static minimization on the resulting averaged minimum-energy structure. 
Here, 'static' means running energy minimization on the electron distribution without chang- 
ing the positions of the nuclei. From this run, we obtained the total energy of the system, 
-Esyst- We then took the same structure and removed the adsorbate to obtain the clean-slab 
structure. Again, we ran selective dynamics on this slab structure to obtain -Egiab- To get 
the energy of an isolated halide atom, we also ran static minimization on an isolated halide 
atom to obtain -Ehai- We define the adsorption energy per supercell per site as 



ads 



-Esyst — -E'slab ~ 2A^i?hal 

2iV;Z ■ 



(2) 



Here, A^^site is the number of sites on one surface of the metal slab, and is the number 
of halides on each side of the slab. In Fig. |2] we show i^ads as a function of Q for both 
systems. We emphasize that i^ads contains the lateral interaction energy and is different 
from the single-particle binding energy E\^. The relation between these two quantities is 
given explicitly in Eq. (IT^ below. 

To understand the surface polarization we need to study the charge-transfer behavior. 
We define the negative of the electron densities from the DFT output as the charge density 
distributions p(x), and we introduce the charge transfer function per adsorbed atom, which 
is defined as follows 22 1 



Ap(f) 



TV 

P(a^) halide- Ag(lOO) ~ p(3;)halide ~ P(a^)Ag(100) 



(3) 
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FIG. 2: Adsorption energy vs coverage. Despite the appearance on this scale, the hnes are, in fact, 
shghtly convex due to the repulsive lateral interactions. The "error bars" in this and subsequent 
figures do not represent statistical errors, but rather estimates of the accuracy of the results, based 
on the convergence studies discussed in Appendix A. 



where p(x)haiide-Ag(ioo) is the full charge density of the adlayer system with adsorbed Br 
or CI on each side of the slab, and p(x)haiide is the full charge density of the pair of isolated 
halide atoms at the same positions as in the halide-Ag bonded system, and p(x)Ag(ioo) is the 
charge density of the Ag(lOO) slab with all atoms at the same positions as in the halide-Ag 
bonded system 23] . After integrating over x and y, this yields the charge transfer function 
per pair of adsorbed atoms, 



Ap{z) = [p(^)halidc-Ag(100) - Np{z\^iide " p{z) Ag{100)] / N . (4) 

From the charge transfer function integrated over the x and y directions, Ap{z), we can 
calculate the surface dipole moment as 

1 r+'' 

p= - y ^ \z\Ap{z)dz. (5) 

Here h = where H is the height of the supercell. The zero point of the coordinate 
is placed at the middle of the supercell. Figure |3] shows the results of the dipole moment 
calculation for Bromine and Chlorine. Here we observe that the magnitude of the dipole 
moment decreases approximately linearly with 6. The surface dipole moment of the energy- 
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FIG. 3: Dipole moment vs coverage. 



minimized clean slab was also calculated and verified to be the same as the surface dipole 
moment of the slab with all atoms at the same positions as in the halide-Ag bonded system, 
thus justifying our procedure. 

Figure m shows the charge transfer function ^p{z) for Br/Ag(100) with 9 = 1/9. In 
this figure positive values indicate electrons being removed, while negative values indicate 
electrons being added. From the figure we see that charge is mostly transferred from the 
surface silver atoms to the adsorbates. Inside the bulk, the charge transfer function indicates 
only minor charge redistribution above and below each of the silver layers. Since the charge 
transfer function is calculated by subtracting the charge distributions of the clean slab and 
isolated adsorbate from that of the adsorbed system, we conclude that this small charge 
redistribution is caused by the adsorption processes. 

Figure |5] shows the charge transfer function per adsorbed atom Ap{z) for Br/Ag(100) 
and Cl/Ag(100) for all coverages. Here we only show half of the supercell since the charge 
transfer function is symmetric in the z direction. Both systems show a similarity in that the 
magnitude and distribution of the charge transfer from the Ag surface are independent of the 
coverage. However, Fig. |5] also shows that while the magnitude of the charge transfer from 
the surface to the adsorbate is independent of the coverage, the resulting charge distribution 
around the adsorbate is not. Indeed, higher coverage results in a more asymmetrical charge 
distribution around the adsorbate. This asymmetry is more pronounced in the Br/Ag(100) 
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FIG. 4: The charge transfer function /S.p{z) for Br/Ag(100) with 9 = 1/9. The vertical short- 
dashed hnes indicate the z-positions of the Ag layers, and the long-dashed ones indicate those of 
the adsorbate ions. 

case, suggesting an important difference between Br/Ag(100) and Cl/Ag(100). Figure [6l 
wfiicli shows the charge transfer function for low and high coverages, illustrates the difference 
more clearly. Here we see that there is no significant difference between Br/Ag(100) and 
Cl/Ag(100) for 6 = 1/9, while for 9 = 1/2 we see a quite significant difference. 

III. DIPOLE-DIPOLE INTERACTION 



In the previous section we have shown that once we have obtained the charge transfer 
function, we can calculate the dipole moment p from Eq. ([5]). Kohn and Lau [24] showed 
that the non-oscillatory part of the dipole-dipole interaction energy between the adatoms 
behaves as 

2PAPB 



V^dip-dip ^^^^^3- 

The novel aspect of this expression is the factor of 2. A qualitative explanation for this 
factor is given in Appendix B. For a more detailed and general treatment we refer the reader 



to Ref. 



24| . With Eq. (|6]), we can calculate 0dip-dip from the surface dipole moment results 
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FIG. 5: Charge transfer functions A/9(z) for all coverages. Only half of the supercell is shown, 
from z = oA to z = 12.235A. Panel (a) is for Br/Ag(100) and (b) for Cl/Ag(100). The dotted 
lines correspond to the z-position of the topmost layer of metal and the dashed lines correspond 
to the z-position of the adsorbates. 



from the DFT as described in Eq. ([5]) as 



2p' 

'dip— dip nnn . j^n 



(7) 



for large R (in our case larger than the nearest-neighbor distance). Here p is the surface 
dipole moment calculated from the charge transfer function (Eq. ([5])), and -Rnnn is the lateral 
distance between a pair of next-nearest neighbor adatoms. 



IV. LATTICE-GAS MODEL 

We use an L X L square array of NsHe = L"^ adsorption sites. Each site corresponds to a 
four-fold hollow site on the Ag(lOO) surface. The energy of this lattice-gas model is 
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FIG. 6: Comparison of charge transfer functions Ap(z) for for Br and CI at two different coverages. 



£■ = - ^ (pijCiCj - 



(8) 



i<j 



Here i and j denote adsorption sites, 0^ is the lateral interaction energy of the pair [ij), 
and i?b is the single-particle binding energy. The sign convention is that (pij < signifies 



repulsive interaction and > favors adsorption 



10| . Sj<j is a sum over all pairs of sites, 



and Nsite is the number of four-fold hollow sites on each side of the slab. For simplicity we 



ignore multiparticle interactions 

Koper ^ has shown that the effects of screening and finite nearest-neighbor repulsion are 
very small. Following his results, we use a lattice-gas model with nearest-neighbor exclusion 
and unscreened dipole-dipole interactions. The distances used in the lattice-gas model are 
Rij = Tija and -Rnnn = v^a, where Rij is the distance between a pair of adsorbates ij, and 
a is the Ag(lOO) lattice spacing. We can then write 



p3 



iV2f 



(9) 
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Thus we have 

7— X] <Pi3CiCj = (f>nnn^9, (10) 

where 



(^2)3 V- CiC 



The adsorption energy defined in Eq. (|2]) is related to the lattice-gas energy of Eq. ([8]) as 

E^ds = -j-^- (12) 

^^site 

This enables us to break down E^ds into its lateral-interaction and single-atom binding parts 
as follows, 

Eads = —(j^nnn^e — E^O, (13) 

where 9 is the coverage (Eq. ([T])) as before. The subscript 9 in signifies that the lateral 
interaction energy is coverage dependent. 

Using the supercell set-up of the DFT, the lateral part of Eq. ([8]) will be the lateral 
interaction energy per supercell surface. We can calculate this energy by extending the 
supercell to infinity in the x and y directions by means of periodic boundary conditions. 
The central supercell is the original supercell, and the image supercells are the supercell 
extensions in the x and y directions. The lateral energy per supercell is the sum of the 
interaction energies of pairs in the central supercell and the lateral energies of pairs of 
adsorbates in the central supercell and adsorbates in the image supercells. Figure [Tl^a) 
shows an example of the lateral energy calculation for = 1/9 for finite A^^site. 

The lateral energy per supercell can be written as 

where C is an arbitrary constant. The above sum can be approximated by the integral 

Se(L) ^ / / —rdrd9, (15) 

which gives us 

S,(L)^^ + C2. (16) 
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FIG. 7: (a)The lateral energy per supercell as a function of the number of adsorption sites (A'site = 
L^). As is increased, the lateral energy approaches an asymptotic value that can be found by 
plotting the energy per supercell as a function of 1/L and extrapolating the graph to 1/L = as 
shown in (b). 



We therefore plot versus 1/L. It is shown in Fig. [Tl^b) that the plot is linear in accordance 
with Eq. (1161) . The correct lateral energy per supercell can then be obtained by fitting 
Eq. flTB]) to the versus 1/L plot and extrapolating to 1/L = 0. The results of this 
calculations for the different coverages are presented in Table HI 
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TABLE I: Lateral interaction energy calculation extrapolated to L — t- oo. is the full lateral 
interaction energy in units of (pnnm the lateral energy between a pair of next-nearest neighbors. 



9 ^e{L^oo) 

1/9 0.10512 

2/9 0.58591 

1/4 0.79822 

1/3 1.53990 

1/2 4.26730 



V. LATTICE-GAS FITTING 

According to our assumption, 0jj is quadratic in p and ~ The ~ l/r"^ part has already 

been calculated in Eg as described in Eqs. ( ITTl) and ( IT4lfT6l) . We also know from the DFT 
results that the dipole moment p is approximately linear in 6 as shown in Fig. [3l Hence, 
based on Eq. ([7]) it is reasonable to assume that we can write 0nnn as {q] 

</)nn„ = A{1 + Bef (17) 

From Eqs. ( fTTl) and ( IT3|) . we have three parameters to be extracted: A, B, and -Eb- In 
Fig. |2] it is shown that i^ads vs 6 is predominantly linear. The linear part is proportional 
to -Eb- The lateral energies contribute to the nonlinear parts which are much weaker, and 
therefore difficult to estimate accurately from a direct three-parameter fit. We therefore 
used the following two-step procedure. As can be seen in fig. [21 the graphs extrapolate to 
Eads{(^ = 0) = 0, consistent with the fact that at a very low coverage the lateral energy 
approaches zero. To obtain the dominant linear coefficient i?b, we first fit a quadratic 
equation to Eadsi^), 

E,ds{0) = ao + a^e + a2e^ (18) 

We extracted the linear part + ai6 and used Oi as our estimate for the linear coefficient 
Eb, finding Eb = 3.059 ± 0.058 eV for Bromine and E^, = 3.371 ± 0.058 eV for Chlorine. 
We then fixed i?b in Eq. ( IT3l) and applied a two-parameter fit to extract A and B, which 
enabled us to calculate 0nnn- (The parameters ao and 02 are complicated functions of A and 
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FIG. 8: The contributions of the lateral interactions to -Eads; shown vs 6. 

B and were discarded in favor of the direct two-parameter fit of the latter.) 

Using E\y from above, we calculate the contribution of the lateral interactions to i^ads as 



c{e) = E^A,{e) + E^e. 



(19) 



In Fig. [8] we plot Eq. f|T9|) . From this figure it is obvious that the lateral energy terms are 
important. Figure [9] shows the fitting results for 0nnn- It is shown in the figure that for Br the 
lattice-gas model obtained by fitting to the adsorption energies from the DFT calculation 
is consistent with long-range dipole-dipole lateral interactions using the dipole moments 
calculated from the DFT charge distribution. This indicates that long-range dipole-dipole 
interactions are dominant in this system. For CI the figures show that the long-range dipole- 
dipole interactions are important but not dominant. 

We further note that for low coverages our estimates of 0nnn for Br are in excellent 
agreement with those obtained by fitting Monte Carlo simulation results for the lattice-gas 
model to electrochemical adsorption isotherms in Ref. js]. However, the DFT results show 
a stronger coverage dependence than obtained from the experimental Monte Carlo fits. The 
experimental fitting results for CI from Ref. js] lie between the two DFT estimates, and all 
three results show approximately the same coverage dependence. 



14 



> 



> 





-0.005 
-0.01 

-0.015 
-0.02 

-0.025 
-0.03 

-0.035 



-0.005 
-0.01 

-0.015 
-0.02 

-0.025 
-0.03 

-0.035 



(a) 



0.1 



(b) 



C3 



0.1 



1 ' r 

Bromine 




G-O Dipole-dipole DFT 

O-^DFTE^^^fit 

□ □ Simulated isotherms fit 



0.2 



0.3 

e 



0.4 



0.5 



1 ' r 

Chlorine 

^ 

[] 

















G-O Dipole-dipole DFT 

O-^ DFT E^^^ fit 

□ □ Simulated isotherms fit 



0.2 



0.3 

e 



0.4 



0.5 



FIG. 9: Comparison of the estimates of (f)nnn from lattice-gas fit to the adsorption energies (dia- 
monds) and from the dipole-dipole interactions (circles). Also shown are the results of fitting room- 
temperature Monte Carlo simulations to electrochemical adsorption isotherms from Refs. 8|, |9| 
(squares), (a): Bromine, (b): Chlorine. 

VI. DISCUSSION 



The lattice-gas model in our study consists of two terms, the lateral interaction term and 
the single-atom binding-energy term. By fitting the lattice-gas model to adsorption energies 
obtained from DFT calculations, we have calculated the total lateral energy of the systems. 
From the charge distribution results from DFT, we have calculated the long-range dipole- 
dipole interaction contribution to the lateral energy terms that falls off as ~ With 
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this assumption, we calculated dipole-dipole lateral interactions by Eq. (Q. 

Apart from the difference of magnitude of the dipole moments between Br/Ag(100) and 
Cl/Ag(100), we find that there are differences in the charge distribution around the ad- 
sorbates between Bromine and Chlorine. This is an indication that there are important 
differences between Br/Ag(100) and Cl/Ag(100). 

For Bromine, we showed that the lateral energy calculations from the DFT charge dis- 
tributions are consistent with the results from fitting the lattice-gas model to the DFT 
adsorption energies. This shows that in the case of Bromine the lateral energy terms are 
dominated by long range dipole-dipole interactions. In the case of Chlorine, the lateral en- 
ergy results from the charge distributions are greater in magnitude than those of Bromine, 
showing that the long-range dipole-dipole interaction in Cl/Ag(100) is important. However, 
in the case of Chlorine, we see less consistency between the two methods of calculations. 
This indicates the presence of significant short-range interactions. 

Our calculations were done in vacuum. We note, however, the overall consistency of 
the vacuum DFT calculations presented here with previous fits of lattice-gas Monte Carlo 
simulations to electrochemical adsorption isotherms. This suggests that our calculations 
might be useful to understand these experimental results, in which water is present, as well. 
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Appendix A: Convergence Checks 

The number of metal layers in our DFT simulation was determined by convergence checks. 
We calculated i^ads for = 1/9 and 1/2, for 5, 7, and 9 layers with exactly the same 
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simulation parameter set-up (energy cutoff, /c-points, the thiclcness of the vacuum regions, 
etc.). From Table |T] we see that increasing the number of metal layers from 5 to 7 changed 
i^ads for Bromine by less than 1 meV for ^ = 1/9 and less than 10 meV for 9 = 1/2. Similar 
observations are also shown in Table [ITl for Chlorine. Increasing the number of metal layers 
from 5 to 7, changed i^ads for Chlorine by less than 2 meV for 9 = 1/'^ and less than 10 meV 
for 9 = 1/2. 

We also calculated the surface dipole moments for 9 = 1/9, and 1/2, for 5, 7, and 9 
layers from the above simulations. The convergence check for dipole moments as shown in 
Table IIIII and IIVI shows that increasing the number of layers from 7 to 9 did not change the 
dipole moment significantly. 

We calculated the percent errors, defined as follows 



PE, 



PEp 



^ads(j) 



X 100%, (Al) 



X 100%. (A2) 

P{3) 

Here, PEe is the percent error for adsorption energies i^ads and PEp is the percent error 
for surface dipole moments. In our calculation j = 5 represents the slab with 5 layers, and 
i = 7,9 represent the slabs with 7 and 9 layers, respectively. These percent errors are also 
shown in Tables Hl UVl 

From these two convergence checks (-Eads and p) we concluded that we need at the very 
least 5 layers of metal, and we decided to use 7 layers. Taking the highest value of -Eads(0 ~ 
Eads{i — 2) from 5 to 7 layers, which is 7 meV, we estimate the error bars for i?ads to be 
AEads = ±10 meV and for p to be Ap = iO.OleA. Error-bar estimates for (pnnn based 
on Ap were then calculated by direct error propagation. Error-bar estimates for E'ads were 
obtained as those leading to a 10% increase in the of the two-parameter fit. 



Appendix B: The Factor 2 in Eq. ^ 



Following Ref. 2J], a qualitative explanation for the factor 2 in Eq. (jH]) can be obtained 
as follows. Consider an adatom A with induced charge qa, at a distance za above the 
plane surface of a semi-infinite conducting medium, located at 2; = 0. The charge-transfer 
function, integrated over x and y, is 
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TABLE I: Convergence check for the Bromine adsorption energy (in units of eV) with respect to 
the number of metal layers. 



BROMINE 



Metal Layers 


Coverage 


Eads 


EadsW - E^s{i - 2) 


PEe 


5 


1/9 


-0.334187984 






7 


1/9 


-0.333754808 


0.0004331 


0.13 


9 


1/9 


-0.339772195 


-0.00601 


1.67 


5 


1/2 


-1.475041628 






7 


1/2 


-1.479772329 


0.00473 


0.32 


9 


1/2 


-1.489228380 


-0.00946 


0.96 


TABLE IL Convergence check for the Chlorine adsorption 


energy (in units of eV) with respect to 


the number of metal layers. 








CHLORINE 










Layers 


Coverage 


Eads 


Eadsii) - Es^{i - 2) 


PEe 


5 


1/9 


-0.371225625 






7 


1/9 


-0.370121449 


0.001104 


0.29 


9 


1/9 


-0.376717001 


-0.005491 


1.48 


5 


1/2 


-1.642027259 






7 


1/2 


-1.649943352 


-0.007916 


0.48 


9 


1/2 


-1.730707884 


-0.080765 


5.40 



^Pa{z) = -QA^iz) + qAS{z - za), (B1) 
where 5{z) is the Dirac delta function. This yields the dipole moment, 

/oo 
zAp{z)dz = -qa- + qaZa = QaZa (B2) 
■oo 

This is the physical dipole created by adatom A. However, the electrostatic potential at a 
point zb, a lateral distance R » za from A, is that of the dipole formed by qa and its 
image charge —qa at z — —za, 
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TABLE III: Convergence check for the Bromine surface dipole moment (in units of eA) with respect 
to the number of metal layers. 



BROMINE 



Layers 


Coverage 


P 


p{i) -p{i - 2) 


PEp 


5 


1/9 


-0.241532 






7 


1/9 


-0.166162 


0.075369 


31.20 


9 


1/9 


-0.178799 


-0.012637 


25.97 


5 


1/2 


-0.124252 






7 


1/2 


-0.118376 


0.005876 


4.72 


9 


1/2 


-0.120788 


-0.002412 


2.78 


TABLE IV: Convergence check for the Chlorine surface dipole moment (in units of eA) with respect 


to the number 


of metal layers 








CHLORINE 










Layers 


Coverage 


P 


p{i) -p{i - 2) 


PEp 


5 


1/9 


-0.267334 






7 


1/9 


-0.199436 


0.067898 


25.39 


9 


1/9 


-0.209752 


-0.010316 


21.54 


5 


1/2 


-0.148922 






7 


1/2 


-0.155043 


-0.006121 


4.11 


9 


1/2 


-0.151803 


0.003239 


1.93 



UA{zB,R) = 2zAqAj:^, (B3) 

for zb > 0. This is equivalent to the potential of a fictitious dipole of magnitude 2zAqA = 
2pA, twice the magnitude of the physical dipole in Eq. flB2p . 

An adatom B with induced charge qs at zb corresponds to the charge transfer function 

ApBiz) = -qB6iz) + qB6iz- Zb), (B4) 

which gives pb = qs^B- 
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The potential energy of the pair of adatoms is then 



/ + 00 
UA{z,R)ApB{z)dz 
-oo 



-qBUA{0,R) + qBUA{zB,R) 

^ZaQaZbQB 



'^^^^ (B5) 



which is Eq. ([6]). In Ref. 2J| it is shown that this result holds in more general situations as 
well, such as jellium and crystalline metals. 
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